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ABSTRACT 


A  laser  reflection  technique  Is  used  to  Investigate  the  relaxation  of 
SlOt  film  stress  which  occurs  during  the  dry  thermal  oxidation  of  SI 
between  700^  and  1000**C.  Included  Is  a  determination  of  the  stress 
distribution  In  the  oxide  by  two  Independent  methods:  1)  measurement  on 
oxides  of  various  thicknesses  from  100  to  800  A,  and  2)  repeated  stress 
measurements  on  chemically  thinned  S10«  films,  viz.  etch  back  analysis 
toward  the  Interface.  Agreement  Is  found  between  these  experiments. 
These  thin  film  stress  measurements  are  achieved  by  the  use  of  ultra- 
thin  SI  substrates  (75  pm).  Essentially,  an  Increase  In  film  stress 
with  decreasing  film  thickness  Is  observed.  Rapid  stress  relaxation  Is 
observed  for  all  temperatures  studied  and  Is  attributed  to  a  time- 
dependent  oxide  viscosity.  The  Influence  of  these  measured  properties 
on  the  kinetics  of  SI  oxidation  Is  discussed. 


INTRODUCTION 


A  number  of  studies  have  reported  the  observation  of  an  Intrinsic 
film  stress  resulting  from  the  thermal  oxidation  of  SI  (1-4).  The 
origin  of  this  stress  has  been  attributed  to  the  120%  molar  volume 
expansion  vhlch  results  from  the  conversion  of  SI  to  S10« .  A  Maxwell 
viscoelastic  flow  model  (2,3,5)  has  been  proposed  to  explain  the 
occurrence  of  Increasing  stress  at  reduced  oxidation  temperatures 
(<950<’C>,  where  higher  oxide  viscosities  preclude  stress  relaxation 
during  normal  oxidation  times.  Since  the  Intrinsic  stress  exists  at 
the  oxidation  temperature.  It  Is  possible  that  this  stress  Influences 
the  oxidation  kinetics.  However,  there  Is  some  disagreement  on  whether 
this  stress  effect  Is  more  Influential  on  the  Interface  reaction 
(2,6,7)  or  the  dlffuslonal  transport  of  Ot  through  a  growing  oxide  (8- 
10).  Also,  other  properties  such  as  SlOt  film  density  (2,11)  and 
electrical  characteristics  (12-13)  may  be  Influenced  by  the  stress  and 
Its  distribution. 

We  report  new  measurements  on  the  SlOa  film  thickness  dependence 
of  Intrinsic  stress  using  a  laser  reflection  technique  (3).  This 
dependence  Is  determined  by  two  Independent  experiments.  In  one 
experiment  the  oxidized  samples  are  grown  to  various  thicknesses,  while 
the  other  uses  oxides  that  are  chemically  etcheo.  Concordant  results 
are  obtained  and  compared  with  present  stress  models.  The  Implications 
of  these  results  on  related  oxidation  models  are  also  discussed. 


EXPERIMENTAL  PROCEDURES 


A.  Sampl*  Preparation. 

All  experiments  were  performed  using  commercially  available  p-type 
(100)  Si  wafers.  Wafer  characteristics  include  a  substrate  thickness 
of  approximately  75  pm  (3  mil),  a  resistivity  of  2  O-cm,  and  an  initial 
diameter  of  2  inches.  Samples  used  for  measurement  were  cut  from  the 
original  wafers  by  cleaving  along  the  major  axis  (100)  to  obtain  1x4  cm 
strips.  A  pre-oxidation  at  high  temperature  (lOSO^C)  was  performed 
prior  to  cleavage  of  the  original  wafers  in  order  to  prevent  subsequent 
contamination  by  Si  dust  during  cleaving.  After  removing  this  oxide  In 
a  concentrated  HF  solution,  samples  were  cleaned  using  a  slightly 
modified  form  of  the  RCA  method  (14)  prior  to  the  experimental 
oxidation. 

Thermal  oxidation  of  the  Si  substrates  was  performed  at  700,  600, 
and  1000*^0  in  a  double-walled  reslstlvely  heated  fused  silica  tube 
furnace.  Oxidation  temperatures  are  reported  with  less  than  2^C 
variation,  and  0*  with  less  than  O. 5  ppm  hydrocarbons  was  used.  The 
Ht 0  content  of  the  was  measured  during  oxidation  at  the  furnace 
outlet  with  a  hygrometer,  and  found  to  contain  less  than  5  ppm  HaO.  A 
flowing  Na  ambient  was  maintained  between  the  inner  and  outer  walls  of 
the  furnace  tube  and  Inside  the  furnace  when  oxidations  were  not  in 
progress.  Film  thicknesses  were  measured  with  an  accuracy  of  better 
than  2%  using  a  precision  manual  ellipsometer  with  a  He-Ne  laser  (6328 
A)  as  the  source.  As  previously  detailed  (2-4),  the  grown  SiOa  was 
removed  from  the  backside  of  each  sample  prior  to  the  strain 


measurement,  and  after  the  measurement  the  substrates  were  rechecked 
for  residual  curvature  with  the  oxide  completely  removed.  Two 
Independent  experiments  were  performed  to  determine  the  distribution  of 
stress  In  SlOa  films.  First,  a  aeries  of  oxides  from  approximately 
100-800  A  were  grown  to  determine  If  any  measurable  stress  relaxation 
occurs  In  the  Initial  thin  regime  of  SI  oxidation.  Stress  measurements 
were  performed  Independently  at  each  thickness  for  the  three  reported 
temperatures.  Since  the  measured  stress  la  the  average  oxide  stress 
Integrated  over  the  thickness  (IS),  this  experiment  establishes  the 
distribution  of  stress  from  the  Interface  to  the  outer  surface  when  a 
series  of  different  thicknesses  are  measured.  Secondly,  oxides  were 
grown  to  about  800  A  and  chemically  thinned  In  a  buffered  HF  solution 
( NH* F; HF/50 : 1 ) .  Stress  was  measured  after  each  etch  on  the  same 
samples.  The  comparison  of  stresses  from  these  two  techniques  yields 
Information  relative  to  stress  relaxation  during  growth,  since  the  thin 
oxides  obtained  In  the  first  experiment  are  exposed  to  Jess  high 
temperature  than  those  obtained  from  the  etch  back  method. 

B.  Strain  measurements. 

A  previously  reported  (3,4)  parallel  laser  reflection  technique 
was  used  to  study  film  stress  as  a  function  of  oxide  thickness,  and 
only  the  essential  features  are  repeated  here.  This  technique  Involves 
a  measurement  of  the  change  In  the  radius  of  curvature  (R)  of  the  SI 
substrate  Induced  by  stress  from  the  SlOt  film.  The  curvature  change 
Is  obtained  from  measurements  of  the  change  of  the  separation  of  two 
parallel  laser  beams  reflected  from  the  film-substrate  surface.  The 


total  film  stress  can  be  calculated  using  a  form  of  Stoney's  equation 
(16)  : 

o,  =  Et. • /6( l-v)t, R  (1) 

where  t.  and  tr  are  the  substrate  and  film  thicknesses,  respectively, 
and  E  and  v  are  Young's  modulus  and  Poisson's  ratio,  respectively,  and 
are  reported  elsewhere  (4).  Since  at  Is  measured  at  room  temperature. 
It  contains  a  thermal  expansion  component  Ot n  In  addition  to  the 
residual  Intrinsic  stress  o,  .  at  n  Is  the  result  of  a  thermal  expansion 
mismatch  between  SI  and  SIO* ,  which  upon  cooling  from  oxidation  to  room 
temperature  yields  a  stress  that  may  be  calculated  (17)  and  subtracted 
from  a,  to  obtain  at  .  The  new  stress  results  reported  herein  Include 
the  extension  of  past  measurements  0>  to  S10«  film  thickness  below 
~a00  A  in  order  to  determine  the  stress  distribution  in  the  oxide. 

This  extension  Is  accomplished  through  the  use  of  thin  SI  substrates 
( A.7S  pm),  so  as  to  yield  measurable  R  values  for  small  thicknesses. 

The  reported  data  In  both  experiments  was  usually  the  average  of  at 
least  3  samples,  each  with  multiple  measurements,  and  resulted  In  an 
average  error  of  -vO.  SxlO*  dyn/cm*  .  The  etch  back  results  showed  a 
larger  scatter  from  the  curves,  which  may  be  attributed  to  the  repeated 
handling  of  the  same  samples  during  the  series  of  etches  and 
measurements. 

RESULTS  AND  DISCUSSION 
A.  Viscous  flow. 

Figures  1  and  2  show  the  dependence  of  Intrinsic  film  stress  on 
oxide  thickness  for  three  temperatures  and  a  range  of  oxide  thicknesses 
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from  100-800  A.  The  data  In  Fig.  1  was  obtained  by  growing  oxides  to 
various  thicknesses  and  In  Fig.  2  from  chemical  thinning  or  etch  back 
of  SlOa .  Both  techniques  yield  essentially  the  same  result  within 
experimental  error.  Considering  that  the  thin  SlOa  samples  used  were 
exposed  to  decidedly  different  thermal  budgets,  the  stress  relaxation 
must  be  fast  relative  to  oxidation  times.  This  Is  a  fundamental  point 
and  will  be  discussed  below.  Also,  there  Is  clearly  a  trend  toward 
larger  stress  for  smaller  SlOa  thicknesses.  This  variation  with 
thickness  Is  smaller  for  low  temperature  grown  SlOa  films  which  exhibit 
larger  stresses.  While  the  data  does  not  extend  to  films  less  than  100 
A,  a  common  extrapolated  value  of  '«.4.  5x10*  dyn/cm*  Is  obtained  for  the 
maximum  Intrinsic  stress  at  zero  SlOa  thickness.  This  finding  of  a 
single  value  for  the  maximum  stress  Is  expected  since  the  origin  of 
stress  Is  due  to  the  constant  molar  volume  change  strain  that  occurs  at 
the  SI  surface  during  oxidation.  The  value  for  this  maximum  stress  is, 
however,  smaller  than  expected  (2)  and  will  also  be  discussed  further 
below.  Finally,  the  stress  values  for  thicker  films  are  near  the 
previously  reported  results  for  thick  SlOa  samples  prepared  on  thicker 
(200  pm)  SI  substrates  (4). 

We  now  return  to  the  Issue  of  stress  relaxation.  According  to  a 
simple  Maxwell  model  for  a  solid,  viz.  a  spring  and  dashpot  In  series, 
we  expect  the  relaxation  of  stress  to  occur  by  viscous  motion  of  SlOa 
away  from  the  SI  surface.  The  strain  rate,  de/dt,  Is  expressed  as: 

de/dt  =  Cl  /r\  ( 2 ) 

where  Is  the  oxide  viscosity  (2).  Hence  for  smaller  viscosities  the 
strain  rate  is  higher  and  relaxation  occurs  more  rapidly.  The 
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equilibrium  viscosity  Is  an  exponent,  jtl  function  of  temperature  as 
shown  by  the  Andrade  equation  (18): 

r\  =  r\»  exp  (En/RT)  (3) 

where  r\a  Is  a  pre-exponentlal  constant  which  represents  the  minimum 
viscosity  and  Er\  the  energy  of  activation  for  viscous  flow.  The 
effects  of  equilibrium  viscosities  will  be  discussed  later.  The 
relaxation  of  stress  due  to  viscous  flow  In  a  Maxwell  solid  can  be 
summarized  by  the  following  equation: 

CTi  (t)  =  o«  exp<-t/T)  (4) 

where  Is  a  pre-exponentlal  constant  (maximum  oxide  stress  directly 
at  the  Interface),  t  Is  the  oxidation  or  annealing  time  (time  exposed 
to  a  given  temperature),  and  x  the  relaxation  time.  Relaxation  times 
have  been  predicted  from  viscosity  measurements  based  on  the  following 
equation : 

T  =  n/G  <5> 

where  G  Is  the  elastic  shear  modulus.  Since  G  Is  a  moderately  weak 
function  of  temperature  (19)  and  r\  a  strong  function  of  temperature, 
the  temperature  dependence  of  the  relaxation  time  Is  strongly  dependent 
on  the  oxide  viscosity. 

From  literature  values  for  G  (20)  and  extrapolated  values  for  r\ 
(20),  relaxation  times  are  estimated  from  equation  (5)  and  presented 
In  Table  1.  It  should  be  noted  that  these  constants  are  obtained  from 
measurements  on  bulk  glass  materials,  and  may  be  In  error  for  thin  S10« 
films.  The  calculated  x's  suggest  that  the  stress  relaxation  Is  slow 
relative  to  normal  oxidation  times  for  low  temperatures.  This  Is 
partly  substantiated  by  the  higher  stress  and  shallower  slope  of  the 


relaxation  curves  for  low  temperature  grown  films.  However,  if  this 
were  entirely  true  then  we  would  expect  the  curves  in  Fig.  1  to  show 
virtually  a  zero  slope  over  the  oxidation  time  periods  investigated. 

In  addition,  the  etch  back  stress  results  of  Fig.  2  would  yield  lower 
stress  values  than  the  individually  grown  samples  of  Fig.  1  in  the  thin 
film  regime.  This  is  due  to  the  fact  that  the  oxidation  rates  are 
greater  in  the  initial  regime  where  thin  films  are  grown  and  thus  less 
exposure  to  high  temperature  is  required  for  the  individually  grown 
thin  oxides.  For  example,  it  takes  7B  min.  to  grow  a  100  A  film  at 
a00°C,  but  to  add  100  A  to  an  already  grown  700  A  film  requires  ~  460 
min  (21).  This  factor  of  6  difference  in  time  should  enable  greater 
relaxation  of  stress  for  the  etched  back  100  A  sample  which  grew  at  a 
substantially  slower  rate.  Also,  it  should  be  remembered  that  the  SiOt 
film  grows  via  the  inward  migration  of  Ot ,  hence  the  newest  oxide  is 
added  at  the  Si-S10«  interface.  Since  the  stress  values  in  Fig.  2  are 
not  lower  than  those  in  Fig.  1,  the  growth  rate  must  have  a  marked 
influence  on  the  relaxation  during  growth.  Given  the  large 
experimental  error  associated  with  the  etch-back  results,  it  is  not  yet 
clear  how  pronounced  this  effect  really  is. 

It  appears  that  the  initial  stress  relaxation  observed  for  the 
grown  SiOi  films  is  faster  than  one  would  expect  from  the  calculated 
t's  reported  in  Table  1,  with  the  result  that  there  is  a  rapid  lowering 
of  the  stress  to  some  initial  value  from  which  there  is  a  slower 
relaxation.  Thus,  either  our  experimental  results  cannot  be  understood 
by  a  simple  Maxwell  model,  since  a  fast  relaxation  is  observed  for  SiO* 
as  it  is  formed  under  potentially  higher  stress  in  the  early  stages  of 


oxidation,  or  a  decidedly  smaller  oxide  viscosity  than  used  In  Table  1 
exists  In  thin  films.  The  fast  Initial  relaxation  may  be  the  result  of 
an  out  of  equilibrium  viscosity  for  the  newly  formed  layer  near  the  Sl- 
510a  Interface.  This  nonequilibrium  situation  cannot  result  from 
temperature  differences,  since  the  oxidation  Is  Isothermal,  but  Is 
likely  derived  from  time-dependent  plasticity  as  a  result  of  changing 
film  thickness.  As  the  growing  layer  moves  away  from  the  Interface, 
the  SlOa  network  relaxes,  and  an  equilibrium  value  for  viscosity  Is 
obtained.  Thus  two  parallel  relaxation  phenomena  are  envisioned:  first 
the  Initially  fast,  temperature  Independent  relaxation  from  a  maximum 
Intrinsic  stress  value  to  a  much  lover  value;  second  the  Maxwell-llke 
relaxation  of  the  well  formed  and  annealed  SlOt  network. 

In  order  to  model  this  situation,  we  can  reconsider  equation  (2), 
which  expressed  the  viscous  flow  rate,  and  provide  a  more  generalized 
viscoelastic  expression  (22): 

de/dt  =  Ao"  (6) 

where  n  is  a  constant.  For  linear  theory  of  viscoelasticity,  n=l, 

A=l/r\,  and  equation  (2)  gives  the  proper  form  for  the  deformation  rate 
(a/r\).  For  1  <  n  <  10,  non-linear  viscoelastic  behavior  results  (22), 
and  equation  (2)  must  be  modified.  Tiller  (6)  has  predicted  a  value  of 
2  for  n  based  on  consideration  of  the  distance  at  which  stress  relaxes 
an  order  of  magnitude,  yielding  a  nonequilibrium  structured  layer  near 
the  Interface.  Other  studies  on  stress  relaxation  In  glasses  have  also 
shown  n=2  (23),  where  greater  stresses  result  In  non-linear  responses. 
An  approach  by  Tan  and  Goesele  (24)  compares  an  Initial  time-dependent 


viscosity  r\i  (t)  In  a  freshly  formed  material  to  a  well-annealed  time- 
independent  viscosity  revising  the  viscous  flow  rate  equation: 

dc/dt  oe  CT/r\»  (7) 

where  r\t  Is  an  average  viscosity.  This  suggests  that  during  oxidation, 
an  average  viscosity  of  the  oxide  over  the  entire  thickness  range  must 
be  considered  rather  than  one  of  true  equilibrium.  A  number  of  authors 
have  suggested  that  the  freshly  grown  oxide  near  the  Interface  Is 
probably  more  fluid  than  the  oxide  near  the  outer  surface  (5,15,24), 
resulting  In  a  nonequilibrium  layer  near  the  Interface.  Also,  the 
growing  film  may  Impose  an  additional  driving  force  to  the  temperature 
activated  viscous  flow,  resulting  in  a  lower  observed  activation  energy 
for  this  process.  This  force  depends  on  the  changing  film  thickness 
and  growth  rate,  viz.  a  decrease  in  film  growth  velocity  as  the  oxide 
moves  away  from  the  surface.  These  Ideas  seem  to  be  consistent  with 
our  experimental  results. 

B.  Kinetic  models  related  to  stress. 

As  stated  earlier,  there  Is  some  disagreement  In  the  literature  on 
whether  stress  primarily  effects  the  Interface  reaction  (2,6,7)  for  SI 
oxidation  or  the  diffusion  of  0«  through  a  strained  oxide  (S-10).  Non- 
Arrhenius  behavior  has  been  observed  for  both  the  parabolic  (kp)  and 
linear  (ki )  rate  constants  (25-27)  that  are  derived  from  the  Deal-Grove 
oxidation  model.  Stress  Is  thought  to  affect  these  constants,  since  a 
break  In  their  activation  energies  occurs  at  approximately  the  viscous 
flow  point  of  the  oxide  (950®C),  where  stress  Increases  for  lower 
oxidation  temperatures. 


The  stress-diffusion  models  (8-10)  consider  a  change  In  the 


diffusion  mechanism  due  to  stress.  Fargelx  and  Ghlbaudo  (8)  have 
considered  the  effects  of  Ct  on  k,  through  a  modified  expression  for 
the  diffusion  (D)  of  oxidant: 

D  =  Do  exp  (-a^V/kT)  <8) 

where  the  pre-exponentlal  constant  Is  the  dlffuslvlty  through  a  stress 
free  oxide,  a  Is  the  level  of  stress  In  the  oxide,  and  Is  the 
diffusion  volume  change  due  to  or.  The  authors  deduce  a  maximum  oxide 
stress  of  4. SxlO' ‘  dyn/cm*  to  account  for  the  difference  In  the 
activation  energies  for  k^  at  low  (<950®C,  2.2  eV)  and  high  <>950®C, 

1.3  eV )  temperatures,  and  propose  It  to  be  constant  over  a  wide 
temperature  range  ( 780-930<*C ) .  We  find  a  maximum  oxide  stress  at  zero 
thickness  an  order  of  magnitude  lower  ('^4.  5x10*  dyn/cm*  )  that  appears 
to  be  Independent  of  oxidation  temperature.  In  agreement  with  the 
origin  of  stress  as  the  constant  molar  volume  change  that  occurs  during 
oxidation.  It  thus  appears  that  the  observed  stress  Is  too  low  and 
therefore  Inconsistent  with  this  model.  However,  the  Influence  of 
stress  on  diffusion  cannot  be  totally  refuted.  Recent  studies  have 
shown  an  Increase  In  oxidation  rate  upon  annealing  thick  oxides  (28), 
presumably  by  relieving  the  compressive  stress  and  Increasing  D.  Also, 
a  crossover  Is  observed  In  oxidation  rates  with  oxide  thickness  for 
different  surface  orientations  (29),  which  appears  to  be  the  result  of 
a  change  In  oxidant  transport  due  to  stress.  At  the  present  time,  none 
of  these  studies  can  be  considered  conclusive,  hence  the  role  of  stress 
on  transport  Is  still  In  doubt,  and  the  presently  proposed  models  are 
not  fully  In  accord  with  the  newest  stress  data. 
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A  different  approach  has  been  taken  by  Irene  (7)  In  describing  the 
effects  of  stress  on  oxidation  kinetics  using  a  revision  of  the  linear 
rate  constant  obtained  from  the  Deal-Grove  oxidation  model  (25).  The 
following  expression  has  been  proposed: 

ki  a  Cs  1  Cl  /rv  ( 9 ) 

where  Ca i  Is  the  areal  density  of  SI  atoms  for  different  SI 
orientations.  In  a  recent  study  <4)  we  showed  the  effect  of  SI 
orientation  on  film  stress.  The  lower  than  anticipated  stress  found 
for  the  (111)  SI  orientation  and  the  high  rate  of  oxidation  suggests 
that  the  above  expression  for  ki  Is  not  correct  as  written.  Relation 
(9)  Incorporates  the  viscous  flow  rate  as  a  factor  affecting  ki .  This 
model  considers  the  Interface  reaction  to  be  controlled  to  some  degree 
by  the  viscoelastic  properties  of  the  oxide  which  relieve  the  volume 
requirements  for  the  conversion  of  Si  to  SIO* .  As  a  result,  a  tensile 
stress  exists  In  the  underlying  Si,  which  may  Increase  the  reactivity 
of  the  SI.  Recently  It  was  reported  (30)  that  the  application  of  a 
tensile  load  to  a  SI  surface  during  oxidation  resulted  In  a  small 
Increase  In  oxidation  rates  for  thin  films  (<200  A).  Our  present 
experimental  data  shows  that  much  of  the  stress  Is  relaxed  to  some 
nearly  constant  value  within  the  first  1000  A  of  oxide  thickness.  It 
has  been  shown  that  k,  may  be  dominant  even  beyond  1000  A  (31),  and 
thus  If  stress  does  Influence  the  oxidation  rate  significantly  It  Is 
likely  to  alter  the  SI  surface  reaction  rate  constant.  Since  we  have 
shown  that  a  high  stress  exists  In  the  early  stages  of  oxidation  even 
at  lOOO^C,  our  results  demonstrate  that  It  Is  possible  for  stress  to 
Influence  oxidation  kinetics  even  at  higher  temperatures  (>950^0. 
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SUMMARY  AND  CONCLUSIONS 

We  have  reported  the  results  of  stress  relaxation  of  dry  thermal 
SlOt  films  on  very  thin  SI  substrates  which  are  used  to  enhance  the 
sensitivity  of  the  measurement.  Our  results  Indicate  a  much  faster 
stress  relaxation  during  oxidation  In  comparison  with  bulk  vitreous 
silica  measurements,  which  may  be  the  result  of  a  reduced  oxide 
viscosity  and  a  resulting  lower  activation  energy  for  viscous  flow.  It 
Is  believed  that  freshly  formed  oxide  near  the  Si-S10«  interface  Is 
more  fluid  than  that  near  the  outer  surface,  and  non-linear  theory  of 
viscoelasticity  must  be  used  to  explain  these  results.  We  propose  that 
a  changing  film  growth  rate  alters  the  viscoelastic  properties  of  a 
growing  oxide,  providing  an  additional  driving  force  to  viscous  flow. 
Since  this  relaxation  occurs  mostly  In  the  Initially  fast  rate  regime 
of  oxidation.  It  Is  believed  that  Intrinsic  film  stress  may  have  a 
stronger  Influence  on  the  Interface  reaction  than  on  the  transport  of 
oxidant  In  the  film,  and  may  be  Important  even  at  temperatures  above 
the  viscous  flow  point. 
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FIGURE  AND  TABLE  CAPTIONS 


Figure  1.  Intrinsic  film  stress  distribution  as  a  function  of  oxide 
thickness  for  Independent  samples  grown  to  various 
thicknesses. 

Figure  2.  Intrinsic  film  stress  vs.  oxide  thickness  for  chemically 
thinned  oxides,  viz.  etch  back  analysis. 

Table  1.  Relaxation  time  constants  predicted  from  viscosity 
measurements. 
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